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ABSTRACT 

A simulation package has been developed for predicting the influence of immersion, i.e. the presence of a uniform 
liquid layer between the last objective lens and the photoresist, on optical projection Uthography. This technology 
has engendered considerable interest in the microlithography community during the past year, as it enables the 
real part of the index of refraction in the image space, and thus the numerical apertmre of the projection system, 
to be greater than unity. The simulation program described here involves a Maxwell vector solution approach, 
including polarization effects and arbitrary thin film multilayers. We examine here the improvement in process 
window afforded by immersion under a variety of conditions, including A = 193 nm and 157 nm, annular 
illumination, and the use of alternating phase shift mask technology. Immersion allows printing of dense lines 
and spaces as small as 45 nm with acceptable process window. We also examine the effect of variations in liquid 
index on the process window and conclude that the index of the hquid must be known to and maintained within 
a few parts-per-million. This has important implications for the temperature control required in future liquid 
immersion projection systems. 

1. INTRODUCTION 

Liquid inunersion Uthography has risen to a promising candidate position for the micro/nanolithography tech- 
nology roadmap for critical dimensions (CDs) down to perhaps 45 nm. In the last few years, several articles^*^'^ 
have been written on various aspects of liquid immersion lithography. Clearly, numerous factors need to be 
investigated to evaluate the viability of this potentially lower cost Uthography technology, including optical, 
fluid dynamic, and fluid property variability concerns. A simulation package, as described here, has been de- 
veloped to address many of the optical issues. The package accounts for high NA imaging, via a full Maxwell 
vector solution approach, including polarization effects, as well as accounting for multilayers of thin film media 
of arbitrary complex index of refraction. This simulation capabiUty was used to examine the process window 
for Une and space structures, based on depth of focus (DOF), for a wide range of conditions: namely, 193 and 
157 nm wavelength projection systems, numerical apertures (NA) ranging between 0.9 and 1.3, and for circular 
and annular iUumination schemes, as well as for alternating phase shift mask structures. Liquid immersion is 
predicted to significantly improve printability between 45 nm and 90 nm CDs, with changes in NA required at 
specific CD junctions to remain on the optimal operating curve. Also examined here are the lithographic effects 
of changes in the index of refraction of the Uquid and the constraints these impose on Hquid temperature control. 

The outline of our article is as foUows. Section 2 describes the physical and mathematical model we imple- 
mented to predict the effects of including a uniform Uquid between the last objective lens and the photoresist on 
the wafer. Section 3 turns to using this simulation package to ascertain the effects on printability of dense lines 
and spaces of various dimensions, when Uquid iuMnersion is taken into account. Section 4 considers the impact 
of uniform changes in real part of the index of the Uquid, such as those arising from a change in temperature. 
Finally, we end this article in section 5 with a few concluding remarks. 




Figure 1. The physical picture of liquid immersion microlithography. The software that has been developed to simulate 
this arrangement can handle an arbitrary number of layers of films in the "thin-film stack" , with each film described by 
a single complex index of refraction. 



2. MODEL OF LIQUID IMMERSION LITHOGRAPHY 

The model of liquid inmiersion lithography we made use of is shown in Fig. 1. The space between the final 
optical element of the projection lens (represented in Fig. 1 by a hemispherical lens) and the resist coated wafer 
is assumed to be filled with a high-index liquid. This eflfect permits exposure with incident light of large obliquity 
that would otherwise be totally internally reflected at the optic-air interface of a dry exposure system.^ A light 
ray of wavelength A in air entering the hemispherical lens of refractive index n\ena at an angle ^lens contributes 
to an image spatial frequency / = niens sin ^lens/'^- For a ray passing next to the rim of the lens pupil, this 
corresponds to a numerical aperture NA = niens sin ^uns = ^Uquidsin^uquid, where nuquid is the Hquid refractive 
index (assumed to be real for the present discussion) and ^uquid is the propagation angle of the ray in the liquid. 
In order for the ray to propagate through the liquid, the numerical aperture must satisfy NA < niens or nuquid, 
whichever is less. 

Suppose the Gaussian focal plane is at a distance D below the optic-liquid interface. The projection lens 
must provide a ray propagating in the liquid at an angle ^uquid (corresponding to image spatial frequency /) 
with a phase shift ^(/), 

^(/) = - ^^^Uquid \/l - Sin^ ^liquid 

= -^nu<,ui<iyC(I^ . (1) 

up to an additive constant, in order for the ray to arrive at the focal plane with the correct phase. If during the 
exposure the hquid refractive index is actually nuquid (as a result of density fluctuation), the phase delay of the 
ray due to propagation from the optic-liquid interface to a plgme a distance D below the interface is <^(/), 



(2) 



This results in a phase error on the focal plane (where D — D) equal to 4>{f) + 
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Thus, fluctuation in the liquid refractive index during exposure will give rise to defocus and spherical aberrations 
of various orders. 

The above model has been implemented in an aerial-image software for Uquid immersion hthography. Vector 
and thin-film interference effects are taken into account^ to allow simulation of the image intensity distribution 
inside the resist layer. High numerical aperture effects are also taken into account by the use of a suitable 
high-NA factor.S'^ 

3. DOF CALCULATIONS 

Our software can calculate the light intensity throughout the entire thickness of the photoresist. However, since 
many resist chemistries are still under consideration for immersion hthography,^ we concentrate on the optical 
effects expected by including a liquid between the last lens surface and the photoresist. Thus, the results reported 
here are for the light intensity just below the surface of the resist, taking into account the influence of the liquid- 
resist interface but not the absorbance or other properties of the bulk of the resist itself. We use the conventional 
aerial image threshold model procedure®'® to estimate the printability of such structures. We anticipate that 
despite the simphcity of this model, the general trends should still hold even when more realistic resist properties 
are later taken into accoimt, 

For the rest of the article, all results will be based on the hght intensity calculated in the top part of the 
photoresist which is taken to be sufficiently thick so as to prevent reflection from the Si substrate. The various 
layers have the following indices: n = 0.6644+2.04i for the underlying silicon, h = 1.52-f 0.03i for the photoresist 
at A = 193 nm, and h = 1.4-|-0.03i for the photoresist at A = 157 nm, and n = 1.47 (water) and n = 1.38+10~®i 
(perfluoropolyether) for the immersion liquids at A = 193 nm and 157 nm, respectively. The immersion Uquid 
thickness is taJcen to be 1 mm at 193 nm and 0.1 mm for 157 nm. These values reflect our best estimates for 
the appropriate parameters given the current state of immersion hthography^ at 193 and 157 nm; they will 
undoubtedly need to be refined as the technologies develop. Still, we anticipate that the general conclusions 
drawn from the present simulations will remain valid. 

Our initial results use contrast calculations to estimate the printabiUty of a given aerial image. Later, we turn 
to exposure-defocus (ED) calculations, following conventions often discussed in B. Lin's works. The contrast, 



was calculated just below the surface of the photoresist. Figure 2 shows curves of contrast versus focus for four 
sets of equal lines and spaces with dimensions of 75, 65, 55, and 45 nm exposed with a binary mask under circular 
illumination with a = 0.6. In each case, curves for NA = 1.1 and 1.3 at 193 nm and NA = 0.9 and 1.3 at 157 nm 
are shown. All exposures except the NA = 0.9 case are made in Uquid immersion. If we take a contrast of 0.3 
as an acceptable level of printabiUty, then a reasonable estimate of the depth of focus (DOF) can be calculated 
from these and other similar plots, as indicated in Fig. 2(c). In Fig. 2(a), all four wavelength/NA combinations 
can print 75 nm L/S, with the largest DOF occurring for the NA = 1.1, A = 193 nm case. As expected, for a 
given wavelength, as NA increases, one can obtain a higher level of contrast when in best focus; however, DOF 
is typically negatively affected. 

A better perspective is shown in Figs. 3(a) and 3(b). At any given feature size, of course one wants the largest 
DOF in order to have as large a process window as possible. Hence, by 'Svalking" along the upper envelope of 
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Figure 2. Contrast vs. focus for 1:1 L/S structures exposed with a binary mask and circular illumination (a = 0.6) 
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Figure 3. DOF vs. feature size with a binary mask and circular illumination for (a) A = 193 nm; (b) A = 157 nm, as 
calculated from the information in Fig. 2. 



each of the curves in Figs. 3(a) and 3(b), we see roughly the desired point at each stage to best take advantage 
of liquid immersion capability. 

Thus, liquid immersion technology offers the ability of increasing the maximum range of NA beyond unity, 
as well as providing a vehicle to select the most appropriate NA value for printing the features of interest, with 
the maximum DOF allowable. In addition, liquid immersion technology can be readily combined with other well 
known resolution enhancement techniques,^^ such as off-axis illumination, optical proximity corrections (OPC), 
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Figure 4. Comparison of contrast vs. focus for annular (ain = 0.4, crout = 0.6) and circular (a = 0.6) illumination, (a) 
65 nm L/S, A = 193 nm; (b) 45 nm L/S, A = 157 nm. 
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Figure 5. DOF vs. focus for annular illumination {ain = 0.4, ^out = 0.6) (a) A = 193 nm; (b) A = 157 nm. 



phase-shift mask (PSM) techniques, and focus-latitude enhancement exposure (FLEX) methods, thereby 
pushing the limits suggested by Fig. 3. 

Figures 4(a) and 4(b) provide an indication of how annular illumination, characterized by inner and outer 
radial parameters of 0.4 and 0.6, respectively,® complements immersion. In Fig. 4(a), we see that this annular 
partial coherence condition improves the DOF for printing 65 nm Hnes and spaces at NA = 1.3, A = 193 nm, by 
about 30%, while at NA = 1.1, the contrast improves to about 30%, at the borderline of printability. Similarly, 
in Fig. 4(b), for A = 157 nm, we see that the consideration of annular illumination, in combination with liqxiid 
immersion and an NA = 1,3, provides borderline printability for lines and spaces as small as 45 nm. Figs. 5(a) 
and 5(b) show the corresponding DOF for the 193 and 157 nm situations, respectively. By comparing with 
Fig. 3, we see that annular illumination provides a good improvement over piu-e circular illumination. 

Similarly, if we lise Uquid immersion in addition to the well known method of alternating PSM, which works 
well for alternating features of spaces and hnes of equal size, we obtain results such as those illustrated in Figs. 6 
and 7. In order to take better advantage of the alternating PSM methods, a smaller partial coherence of 0.3 was 
used in these simulations. In Fig. 6(a), we see that with alternating PSM, a very acceptable printability should 
be obtained for both dry and Uquid immersion at A = 193 nm, with NA as low as 0.9, for printing 65 nm hnes 
and spaces; this represents a considerable improvement over a binary mask with circular illumination (a = 0.6), 
as in Figs. 2(b) and 3(a). Still, of course, liquid immersion provides a larger process window, even for NA=0.9. 
Moreover, with A = 157 nm, as in Fig. 6(b), we see that the use of alternating PSM plus Uquid immersion 
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Figure 6. Comparison of contrast vs. focus for alternating PSM [p = 0.3) and binary mask (a = 0.6). (a) 65 imi L/S, 
A = 193 nm; (b) 45 nm L/S, A = 157 nm. 
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Figure 7. DOF vs. feature size for alternating PSM (a = 0.3) (a) A = 193 nm; (b) A = 157 nm. 



produces a very acceptable DOF for printing lines and spaces as small as 45 nm at A = 157 nm. Finally, Figs. 
7(a) and 7(b) show the DOF for various situations at A == 193 nm and 157 nm, when the alternating PSM is 
used along with these very high NA applications. Tables 1, 2, and 3 sununarize the data contained in Figs. 3, 
6, and 7. 
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Tablel. DOF in nm of 1:1 line/space patterns exposed with a binary mask and circular illumination (cr = 0.6). 
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Table2. DOF in run of 1:1 line/space patterns exposed with a binary mask and annular illumination 

{(Tin = 0.4, aout = 0.6). 
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Tabled. DOF in nm of 1:1 line/space patterns exposed with an altPSM mask and circular illumination 

{a = 0.3). 

The above analysis has also been generalized from simple contrast calculations, to exposure-defocus analysis, 
thereby enabUng DOF predictions to be obtained for isolated structures, as well as more general ID and 2D 
structures. Figure 8 illustrates the general idea. Figure 8(a) contains a plot of exposure versus focus for two 
situations, 75 nm hues and spaces exposed at 193 nm using a binary mask and circular illumination (cr = 0.6) 
at NA=0.9 (dry) and NA=1.1 (Hquid). The DOF improvement and general increase in process window is quite 
evident. For any given exposure latitude, Fig. 8(a) can be used to calculate DOF, as shown in Fig. 8(b). 
Comparing the 75 nm point in Fig. 3(a), obtained via contrast calculations, we see excellent agreement with the 
10% EL point in Fig. 8. Finally, Fig. 8(c) illustrates what happens if we push the technology using alternating 
PSM in addition to liquid immersion. As expected from the earlier contrast calculations, acceptable process 
windows look very feasible at 55 nm dimensions, with NA=1.1 and 1.3 at A = 193 nm. 

4. EFFECT OF LIQUID INDEX VARIATION 

We next turn to a concern regarding Uquid immersion technology that does not exist under normal conditions 
when the index of the vacuiun or gas in the image space is very close to unity. The Uquid may exhibit variations 
in the index as a function of space and time due to a variety of reasons, such as nonuniformity in mixing, 
bubbles, outgassing resulting from the chemical processes occurring in the photoresist, or, our chief concern 
here, nonuniformities in temperature from heating effects. We concentrate on the simplest treatment of such 
concerns, the effect on imaging due to a uniform change in the real part of the index of refraction of the liquid 
throughout the entire liquid. Clearly, other effects such as bubbles or nonuniform heating would cause problems 
not captured by the present simulations; we hope to be able to address such concerns in future work by making 
appropriate extensions to our simiilation capabilities. Nevertheless, the present study at least addresses the 
first order concern that would need to be met to have adequate imaging capabiUty for meeting the printability 
demands expected of the liquid unmersion technique for microlithography. 

Figure 9(a) shows how a uniform change in index throughout the liquid results in each ray being refracted 

differently, thereby resulting in a different image point corresponding to each point in the object space. The 
larger the thickness of the hquid and the larger the change in index of the Uquid, the more important this 
consideration wiU be. 

Figxu*e 9(b) considers the situation of contrast versus focus for 65 nm lines and spaces, with A = 193 nm and 
NA=1.3. Three curves are shown, corresponding to the situations where the change in the real part of the index 
of the Uquid when the index varies imiformly by ±10 ppm from its nominal value. From this, we can calculate 
the corresponding reduction in the DOF. 

Figures 10(a), 10(b), and 10(c) examine the effect of this uniform change in index on the contrast versus focus 
plots previously considered. Figure 11(a) shows DOF versus the change in index for the circular illumination 
case, when A = 193 nm, NA=1.3, and the Uquid thickness is 1 mm. The slopes of the curves for the cases shown, 
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Figure 8. (a) Exposure-Defocus (ED) plot for 75 nm L/S features at A = 193 nm with binary mask and circular 
iUumination {a = 0.6) at NA = 0.9 (dry exposure; top three curves) and NA = 1.1 (immersion; bottom three curves). 
The middle curve in each set of three represents the dose required to obtain the desired line width while the top and 
bottom curves represent the doses for ±10% variations in line width. The process window is considerably larger for 
the NA=1.1 case, (b) DOF vs. exposure latitude (EL) for 76 nm L/S, A = 193 nm, where EL is defined here as the 
exposure difference between the ±10% lines, divided by the exposure at nominal printing, (c) DOF vs. EL for 55 nm 
L/S, A = 193 nm, using an alternating PSM mask and a = 0.3 illumination. 
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Figure 9. (a) A uniform change Sn in the index of refraction of the Uquid leads (to first order) to a defocus AD. (b) 
The effect of 6n on the contrast for 65 nm L/S, at A = 193 nm, NA = 1.3, cr = 0.6. The hatched areas represent the focus 
range over which contrast > 0.3 for each Sn, illustrating the loss of DOF due to uncertainty m Sn, 



of 60, 65, 70, and 75 nm lines and spaces, are fairly similar; of course, the smaller the line/space dimension, the 
lower the DOF of the structure for any given 5n. Figure 11(b) shows our prediction for the DOF variation as 
a function of the change in index, for the 65 nm situation, for all three cases of circular illumination, annular 
illumination, and alternating PSM. Here also, the slopes of the three curves are similar. Again, though, for any 
given Sn, the DOF of the circular illumination is the smallest, then that of the annular, and then that of the 
alternating PSM. 

The main conclusion to be drawn from Fig. 11 is that the liquid index must be matched to the design index 
to ±10 ppm or better. If the mismatch is due only to the temperature control of the liquid, then the latter must 
be within ±0.1 °C (assuming dn/dT « -100 ppm/°C for water at 193 mn).^^ 



5. CONCLUDING REMARKS 

The simulation package produced and described here has enabled us to examine the expected improvements in 
printability by using liquid inmaersion microlithography. Clearly, the advantages to be gained are considerable, 
since it appears that acceptable printability of lines and spaces with dimensions as small as 75 nm are achievable 

with binary masks and circular illumination at A = 193 mn, with DOF ^ 200 nm [Fig. 3(a)]. When combined 
with annular illiunination, dimensions as small as about 68 nm with DOF ^ 200 nm, still at A = 193 nm, appear 
achievable. When liquid immersion is combined with an alternating PSM method, then dimensions less than 
50 nm appear to be achievable with this same DOF [Fig. 7(a)]. Other resolution enhancements techniques, such 
as line biases, attenuated PSM, serifs, etc., should also work well since their eflPects are additive to those of liquid 
inunersion. 

In future work, we anticipate addressing further nonuniformity issues, as well as analyzing 2D imaging 
conditions, such as contact holes, islands, elbows, etc. With our present simulation package, we can analyze 2D 
imaging situations, as well as the problem where the index of the film varies along the optic axis, but is constant 
within any given plane perpendicular to the optic axis. This scenario will arise as the photoresist is heated during 
exposures, creating a temperature gradient in the liquid along the optic axis.^^ The more complicated issues of 
lateral variation of index, and of extreme changes such as those arising from bubbles in the immersion liquid, 
are more difficult to simulate. However, we anticipate attempting to estimate the effects of such variations in 
our future simulation work. 




Figure 10. The effect of 6n on 65 nm L/S patterns for (a) binary mask, circular illumination (cr = 0.6); (b) binary mask 
with annular illumination (ain = 0*4, aout = 0.6); (c) alternating PSM mask with circular illumination (a — 0.3). 
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Figure 11. Sensitivity of DOF to liquid index variation at A = 193 nm, NA = 1.3, and 1 mm liquid thickness, (a) 
Circular illumination (a = 0.6) and binary mask at different L/S dimensions; (b) Different illumination and masks at 
65 nm L/S. 
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